Theory of the different photoemission spectra of metallic and insulating Ceo 

compounds 
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Metallic KsCeo shows pronounced structure and a sharp Fermi edge in integrated photoemission 
spectra (PES) while the insulating K4C60 and KgCeo phases display only a broad structureless 
peak. We find that both types of spectra can be explained by the coupling to the optic vibrations 
of the /C"^Q ionic lattice. This is suppressed in K3C60 due to metallic screening but is strong 
in the insulating phases. We use the non-crossing approximation to calculate the density of states 
(DOS) of electrons in K3C60 coupled to the intramolecular Hg modes in good agreement with the 
experiment. For K4C60 and KeCso strong coupling to the low energy optic Vi'^ /G^q modes controls 
the DOS and yields broad peaks in the PES. A moment expansion is used to calculate the position 
and width of these peaks which agree well with the experiment. 

PACS numbers: 63.20.Kr, 73.61.Wp, 79.60.Fr 



Alkali-metal-doped Cgo has been investigated exten- 
sively by photoemission spectroscopy (PES) in bulk sys- 
tems (Ref. Q, H and references therein). Recently, 
Yang et al. measured the electronic band dispersion 
of K3C60 in a monolayer system on Ag(lll) by angle 
resolved PES (ARPES) 0]. By comparing angle in- 
tegrated spectra (AIPES) from several experiments a 
generic feature emerges which is the striking difference 
between the spectra of metallic K3C60 and the insulating 
K4C60 and KgCgo (the former is a Jahn- Teller distorted 
insulator, the latter a band insulator): KaCgo shows a 
metallic spectra with a sharp Fermi edge and distinct 
structures whereas the insulating phases display a broad 
Gaussian-peak. In the present work we address this dif- 
ference. The key ingredient is the large change in cou- 
pling strength of the low-energy [oj « 10 meV) optic 
vibrations of the K"*" ions. These modes cause a net shift 
of the Ceo molecular levels and their interaction is effi- 
ciently screened in metallic K3C60 where it is reduced 
by a factor 10""* — 10~^ 0,0- However, in the insulat- 
ing phases, where metallic screening is absent, they cou- 
ple strongly {g m 70 meV, see below) due to the direct 
Coulomb interaction of the ionic charge with the photoe- 
mission hole. The resulting physics is best illustrated by 
the toy- model H — lo a^a + g c {a^ + a) where a single 
electron is coupled to a harmonic oscillator. The corre- 
sponding photoemission spectrum is a Poisson distribu- 
tion P(e) = v'^e~^ ln\ 8{t-\-\n-~v\jS) where v = (g/uj)'^ 
is the average number of excited phonons emitted dur- 
ing the photoemission process. In the insulating phases 
J/ is very large {v ~ 50) and as a result the spectrum be- 
comes incoherent resulting in a Gaussian-like shape. In 
the following we treat the metallic and insulating phases 
separately. First, using the non-crossing approximation 
(NCA), it is shown that the spectrum of metallic KaCgo 
is dominated by the coupling to the intramolecular Cgo 
modes. Second, we use moment expansion to calculate 
position and width of the Gaussian-like spectra of K4C60 



and Kg Ceo. 

The only modes that couple strongly (A « 1) in KaCgo 
are the intramolecular Hg modes which cause a splitting 
of the threefold degenerate LUMO's of Ceo (See Ref.'?' for 
more details) . We treat these modes by NCA ^ and ne- 
glect the on-site Coulomb interaction ([/ sa 1 eV in bulk). 
This is justified in a monolayer adsorbed on Ag(lll) 
where U is reduced by the proximity of the metal 
Here we restrict our attention to such systems. A similar 
calculation was performed by Liechtenstein et al. who 
showed that the width of the plasmon in KsCeo can be 
explained by electron-phonon coupling The Hamil- 
tonian which describes the conduction band of KaCgo 
coupled to the Hg modes is (setting h — 1) 
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The first term is the tight-binding band in standard no- 
tation formed by the 3-fold degenerate (ti„) LUMO's of 
CgQ. The sum is over the orbitals n and to, the lattice 
sites j and the nearest neighbors 5. As the spin ori- 
entation is preserved in the Hamiltonian, explicit sums 
over spins are dropped throughout. The second term is 
the energy uo^ of the eight five-fold degenerate Hg vi- 
brational multiplets . The indices v and k denote the 
multiplet and the mode respectively. The last term de- 
scribes the electron-phonon coupling. The phonon en- 
ergies uj^ and the coupling parameters g^ are disper- 
sionless. The parameters g^ are related to the par- 
tial coupling constant by gl = luj^K/N{0) where 
A^(0) is the density of states per spin and molecule 
Values for Wi, and Xi, where taken from Ref. El The 
structure of the coupling is given by the coefficients 

Cnm = \/|(-l)"(2,A:|l,-m;l,n) where (...) is the 
Clebsch-Gordan coefficient pjj. The normalization is 
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such that J2knm (C'nm) = 3. In the NCA, the electron 
self-energy is determined self-consistently by the lowest 
order self-energy diagram with the interacting Green's 
function. In the present problem the Green's function 
Gnm(w, k) and self-energy E„m(w,k) are 3x3 matri- 
ces. The free phonon propagator Dlf.{uj,q) — D1{lj) 
depends only on frequency and the multiplet. As a 
consequence, only the local part of the Green's func- 
tion G^°^ = ^/NJ2k Gmyii'-o, k) enters which renders the 
self-energy E„m(w,k) = S„m(w) local as well. This is 
a consequence of including only non-crossing diagrams. 
The evaluation of the basic diagram can be simplified 
further when the symmetry of the lattice is used. The 
momentum-independent Green's function Gj^^ and self- 
energy have to be invariant under all symmetry 
transformations which belong to the point-group of the 
lattice. In particular, in a cubic environment (such as a 
fee lattice) 



i{uj) = T.{uj)dnrn, (2) 



where G^°''{uj) and are scalars. On the surface or 

in a monolayer the symmetry is lower than cubic and 
the Green's function has additional off-diagonal parts. 
However, the corresponding corrections were found to be 
small (1 % or less) and therefore it is an excellent ap- 
proximation to use the Green's function and self-energy 
as given in 0. This yields the scalar equation 



u). (3) 



Thus, the problem is simplified to a single band interact- 
ing with a discrete set of phonon modes. Relation ||2Jl, 
which is an equation for the self-energy, is solved it- 
eratively using the advanced Green's function. In 
the band-structure only enters via the density of states 
(DOS) and we chose a generic square DOS with a width 
W = 0.5 eV. Using other bare DOS revealed that the 
interacting DOS depends only weakly on the form of the 
bare DOS. The result for half-filling (/i = 0), which cor- 
responds to K3C60, is shown in Fig. The interacting 
DOS shows an overall structure, such as a dip at 0.2 eV 
and a second hump at 0.4 eV, which agrees well with the 
AIPES of the monolayer system Q. 

The insulating compounds K4C60 and KgCgo differ im- 
portantly from K3C60. As discussed in the introduction, 
in the absence of metallic screening, the coupling to the 
K+ modes is strong and the ground state of a hole cre- 
ated by photoemission will be polaronic. This type of 
physics is not captured by NCA. On the other hand, the 
ground states of the insulating phases are rather simple: 
KgCeo is a trivial band insulator, whereas in K4C60 there 
are four localized electrons on each Cfin which form a sin- 
glet due to Jahn- Teller distortions [ij, ll9j . If the ground 
state 1^0) is known, the moments fik = J rfee'^P(e) of 
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TABLE I: Frequencies and coupling constants for the vibra- 
tional modes in K„Ce,o (all energies are in meV). Parameters 
for the intramolecular Hg and Ag modes were taken from 
Ref . 11^ . The coupling constant for the lattice vibrations were 
calculated in this work. In K3C60 the couphng to the Ag 
modes and the vibrations of the ionic lattice is efficiently sup- 
pressed by metallic screening, a) applies for K4C6o. b) apphes 
for Kg Ceo. 
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FIG. 1: Solution of the NCA at half filhng (^ = 0) which 
corresponds to K3C60. Upper panel: Occupied part of the 
interacting DOS. Lower panel: Advanced self-energy. Inset: 
Occupied DOS (solid line) convoluted with a Gaussian (cr — 
10 meV) and compared to the experimental spectrum of the 
monolayer (dashed line) Q- 



the spectrum P{e) (normalized to 1) can exactly be cal- 
culated by evaluating the expectation values 
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The sum is over all sites j and orbitals n and N is the 
total number of the electrons. H is the full Hamiltonian 
and includes all vibrational modes as well as Coulomb in- 
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FIG. 2: Gaussian fit (solid curves) to the photoemission 
spectrum using the first and second moment as calculated 
in (|1J and plotted with respect to the chemical potential. The 
dashed curves are the experimental data (Ref. for KeCeo, 
Ref. for K4C60) 



teractions. Although relation Q is exact, the reconstruc- 
tion of a distribution from a finite number of known mo- 
ments is an ill-defined problem if the overall shape of the 
distribution is unknown. However, as was argued above, 
the strong coupling to the low-energy optic modes yields 
a very incoherent, and therefore Gaussian-like spectrum. 
Hence, for this physical reason the distribution should be 
well approximated by a Gaussian which is determined by 
the first and second moment as given by The proce- 
dure which follows is, first, to determine all contributions 
to H and, second, calculate the moments by 

The additional phonon contributions to H are modes 
which cause a net shift of the molecular orbitals 
and therefore are no longer screened in the insulating 
phases |^ 0. These modes are the intramolecular Ag 
modes and the vibrations of the ionic lattice. In princi- 
ple, lattice vibrations also couple via a change in hopping, 
but the corresponding coupling is much smaller 0^. In 
the case of the two intramolecular Ag, the frequencies 
Wi, and coupling parameters gi, are dispersionless. As for 
the Hg modes, the latter is given by gl = ^u}i,Xi,/N(0) 



(see above) and values were taken from Ref. [Ij. The 
coupling to the ionic lattice has been studied much less 
extensively, mainly because it is negligible in the super- 
conducting KsCgo. Here we consider the coupling due 
to the Coulomb interaction between the ionic charges. 
The mass ratio Mq^^/M-k — 18.4 is big which allows 
to separate the lattice vibrations into optic dispersion- 
less if-modes and acoustic modes. The frequencies of 
these modes were measured by EELS 0| where it was 



observed that K+ ions close to the surface have sub- 
stantially lower frequencies. As photoemission is surface 
sensitive, we use these values which are lok = 8.9 nieV 
and lok = 10.9 meV for K4C60 and KgCgo respectively. 
In what follows, only the averaged coupling constant 
g\ — A^s^^ X]qa ISqaP cuters (the sum is over all optic 
K-modes a) which is given by (with h exceptionally in- 
cluded) 



2Mkujk 
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The sum runs over all K"*" ions r and is the electric 
field at the positions r and caused by an additional hole 
on the Geo molecule at the origin. The sum in |(SJ| de- 
pends on the lattice. KgGgo has a bcc lattice with a cubic 
lattice constant a — 11.39 A where each Geo molecule is 
surrounded by 24 K+ ions located at (0, 0.5, 0.25)a flij . 
Gonsidering the bare coupling of the closest by K"*" ions 
yields gx — 120 meV. Taking into account the pola riz- 
ability of the Ggg^ ions by multipole expansion |17| re- 
duces the electric fields entering (jSJ by 40 % and yields 
a coupling constant gx = 72 meV. In K4G60 distances 
between Cqq and are almost the same as in KgGgo, 
however, every Cqq ion is surrounded by 16 K"^ ions in- 
stead of 24. This reduces gK in K4G60 by a factor \/2/3 
with respect to Kg Ceo- Taking also into account the dif- 
ference in ujk we find gK = 65 meV for K4Ceo- Coupling 
to the acoustic modes is more involved because both ujqa 
and gqa are q dependent. We used a simple spring model 
parametrized by a phonon frequency of 5 meV at the Bril- 
louin zone boundary This yields an averaged cou- 
pling constant ga = 10 meV. In addition we find an aver- 
age frequency uja = Sqa ^qalSqa/ffaP = 3.8 meV. 
Finally, the Coulomb interactions should also be included 
in H. However, they vanish in the case of KeCeo where a 
single hole is created in a full band. In K4Geo there is a 
contribution from the on-site Hunds-rule coupling term 
which involves an exchange energy K — 50 meV (see 
Ref. I23 for a detailed description) . This term leads to an 
additional increase of the width of the K4Ceo spectrum. 

In the following we will discuss the results for KeCeo in 
more detail than those of K4Geo because the ground state 
of the latter involves Jahn- Teller distorted molecules and 
a detailed discussion would go beyond the scope of this 
letter. As mentioned above, the ground state of KeCeo 
is a full band and trivially given by |\E'o) = Hjn ^]n\'^^'^)- 
Note that jv^o) doesn't have any phonon excitations. 
Coulomb terms can be neglected for KeCeo and the 
Hamiltonian H — i/kin + Hp + Hep consists of the kinetic 
energy, the phonon energies and the electron-phonon cou- 
pling terms. As in the previous section we assume a 
quadratic bare DOS po with a width W — 0.5 eV and 
centered around zero. Using relation to evaluate the 
moments of the photoemission spectrum one finds that 
/ii = 0. The second moment is fj,2 = ^2(^0) + J2x9x 
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where ^2(^0) = W/{2\/3) is the second moments of 
the normahzed square DOS pq. The sum is over all 
phonon modes x. Similarly, the third moment is fi^ — 
— ^x'^x9x- Using the parameters as listed in Tab. ^ 
one finds a = ^ = 0.229 eV and ^.^/a^ = -0.319. 
Usually, the photoemission spectrum is plotted with re- 
spect to the chemical potential which, per definition, is 
p. = EoiN) - Eo{N - 1) = -Eo{N - 1) and therefore 
corresponds to the polaron ground state energy. In the 
small polaron limit, where the hole is localized on a sin- 
gle molecule, the chemical potential is given by the re- 
laxation energy of the phonon degrees of freedom which 
couple to the hole: 



9l 



§1 



UK 



0.599 eV, 



(6) 



Eo{C\q) is the Jahn- Teller ground state energy of an elec- 
tron in the LUMO of an isolated Cgo molecule interacting 
with the intramolecular modes. By particle-hole symme- 
try, this is the same as Eq{C\'^) and was calculated nu- 
merically in Ref. The last two terms are the energy 
gain due to the ionic lattice distortions. It must be noted 
that the total contribution of the K-modes to the chem- 
ical potential is proportional to , because g^- itself 
is proportional to (relation ©). Hence, /i depends 
sensitively on ojk- Note that this estimate is a lower 
bound for the chemical potential, because in the insulat- 
ing KgCgo it may lay everywhere in the bandgap rather 
than on the top of the filled band. In Fig. 12 a Gaussian 
of width a and shifted by is plotted. The result com- 
pares well with the experimental curve from bulk mea- 
surements although the width is somewhat too small ■ 
Other fitting functions which accounted correctly for the 
non-zero third moment were tested as well which, how- 
ever, yielded only slightly different curves than the Gaus- 
sian shown here. 

The case of K4C60 is more complicated, due to Jahn- 
Teller distorted molecules where the three-fold degen- 
erate LUMO is split into a singlet and a doublet. In 
K4C60 the doublet is lower in energy and occupied by 
four electrons. Hence, they form a full band in the solid 
which results in the insulating state. The problem of 
the Jahn- Teller distortion in Ceo , which cannot be solved 
analytically, has been extensively studied by various ap- 
proaches In addition, we developed a variational 
wavefunction for Cg,^ which yields a ground state energy 
in agreement with exact diagonalization results [l9|. In 
order to calculate the expectation values in Q we use 
this variational ground state (details of the calculation of 
the will be published elsewhere) . We find a first moment 
/xi = —0.174 eV which is due to the energy gain of the 
Jahn- Teller distortion. The chemical potential is again 



given by ©, except that — _Bo(CgQ ) has to be replaced 
by £;o(Ce^) - Eo{GI^) which yields /i = 0.422 eV 
Finally, the second moment is cr = ^/Ji2 = 0.244 eV 
which is almost the same as in KgCgo- In Fig. |21 the 
corresponding Gaussian is plotted and compared to the 
experimental data from bulk measurements Again, 
good agreement is found. 

In conclusion, we showed that the different photoemis- 
sion spectra in metallic and insulating K„C6o are due to 
a large change in the coupling strength to the low-energy, 
optic K+ modes. Theoretical calculations for both cases 
yield good results. 
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